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Preface 

This  study  was  made  to  determine  quantitative  gamma 
radiation  rates  from  the  products  of  a  specified  fission 
device.  The  .1  sec  to  2  min  time  after  fission  period  was 
of  special  interest.  This  period  is  of  the  most  interest 
in  weapon  fireball  radiation  calculations. 

I  would  like  to  thank  Dr.  Charles  J.  Bridgman  for 
his  guidance  and  interest  in  this  study  and  Mrs.  Audrey 
L.  Crosby  of  the  Vallecitos  Nuclear  Center  for  her  assist¬ 
ance  in  obtaining  needed  data  references. 

Ronald  B.  Drinkwater 
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Abstract 

An  analytical  study  was  made  to  determine  the  gamma 

radiation  produced  by  the  products  of  a  fission  device. 

The  radiation  during  times  from  .10  sec  to  2  min  after 

fission  was  of  special  interest.  Fission  yield  data  were 

235  238 

gathered  for  5  fuels  (fission  spectra  U  ,  U  ,  and 

ooq  235  22ft 

Pu  ;  and  14  MeV  spectra  U  and  U  ) .  Decay  data 

were  gathered  for  the  fission  products  of  these  fuels. 

A  computer  program  was  written  to  process  the  data  which 

produced  quantitative  energy  and  number  rates  of  interest 

in  fission  product  gamma  radiation  analysis.  Results 

could  be  obtained  for  a  device  of  specified  yield  and  fuel 

composition.  The  program  results  were  not  valid  at  times 

shorter  than  1000  sec  after  fission  because  gamma  decay 

information  w  s  not  available  for  many  short  lived  products. 

However  an  approximation  based  on  the  long  time  results  is 

presented  for  use  in  obtaining  short  time  total  gamma  energy 

rates. 
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GAMMA  RADIATION  FROM 
FISSION  PRODUCTS 

I .  Introduction 


Purpose 

This  study  was  made  to  determine  the  gamma  radiation 
given  off  by  the  fission  products  of  a  specified  device  as 
a  function  of  time  after  the  fission  process  occurred.  The 
study  was  primarily  designed  to  produce  results  for  use  in 
the  study  of  delayed  nuclear  weapon  effects.  The  results 
could  be  used  in  reactor  product  radiation  studies,  espe¬ 
cially  for  severe  excursion  analysis. 

Aircraft  flying  in  the  vicinity  of  a  nuclear  weapon 
detonation  may  not  be  harmed  by  prompt  weapon  effects, 
however  their  flight  paths  may  pass  through  the  radio¬ 
active  weapon  products.  The  gamma  radiation  from  these 
products  may  cause  damage  to  the  aircraft  and  its  crew. 
Elementary  buoyant  body  calculations  indicate  a  weapon 
fireball  would  rise  to  the  tropopause  (above  typical  air¬ 
craft  altitudes)  in  approximately  2  min.  If  the  aircraft 
were  to  enter  the  fireball  in  less  than  0.1  sec  after 
weapon  detonation,  prompt  weapon  effects  would  most  surely 
have  destroyed  it.  Thus,  the  time  period  of  0.1  sec  to  2 
min  after  fission  was  of  special  interest. 

Terms  Defined 

Two  terms  used  frequently  in  this  report  must  be 
clearly  defined;  "Fission  Fragment"  and  "Fission  Product". 
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For  the  purposes  of  this  report ,  a  fission  fragment  is  an 
atom  that  is  formed  directly  by  the  fission  process.  A 
fission  product  is  an  atom  that  is  either  formed  by  the 
fission  process  or  present  as  the  result  of  fission  frag¬ 
ment  decays. 

Overview 

Fission  yield  data  and  decay  data  on  fission  products 
were  gathered  and  compiled  for  computer  processing.  A 
computer  program  was  written  to  process  the  data.  The 
program  output  included  both  energy  and  number  rates  (see 
Table  7} .  The  17  gamma  energy  rates  (rates  2-18)  are  the 
same  as  17  of  the  18  gamma  groups  used  in  an  AFIT  prompt 
effects  code  (SV  Code).  Group  I  (8.00-10.00  MeV)  was  not 
included  because  no  gammas  of  this  energy  were  encountered 
in  this  study. 


GNE/PH/74-3 


Table  I 


Program  Output  Rates 


Race  Index 


Definition 


1  Total  Activity  (Decays/sec) 


2 

6.50-8.00 

MeV 

Gamma 

Rate 

(MeV/sec) 

3 

5. 00-6. 50 

MeV 

Gamma 

Rate 

(MeV/sec) 

4 

4.00-5.00 

MeV 

Gamma 

Rate 

(MeV/sec) 

5 

3.00-4.00 

MeV 

Gamma 

Rate 

(MeV/sec) 

6 

2.50-3.00 

MeV 

Gamma 

Rate 

(MeV/sec) 

7 

2.00-2.50 

MeV 

Gamma 

Rate 

(MeV/sec) 

8 

1.66-2.00 

MeV 

Gamma 

Rate 

(MeV/sec) 

9 

1.33-1.66 

MeV 

Gamma 

Rate 

(MeV/sec) 

10 

1.00-1.33 

MeV 

Gamma 

Rate 

(MeV/sec) 

11 

.80-1.00 

MeV 

Gamma 

Rate 

(MeV/sec) 

12 

.60-  .80 

MeV 

Gamma 

Rate 

(MeV/sec) 

13 

.40-  .60 

MeV 

Gamma 

Rate 

(MeV/sec) 

14 

.20-  .40 

MeV 

Gamma 

Rate 

(MeV/sec) 

15 

.20-  ^0 

MeV 

Gamma 

Rate 

(MeV/sec) 

16 

.10-  .20 

MeV 

Gamma 

Rate 

(MeV/sec) 

17 

.05-  .10 

MeV 

Gamma 

Rate 

(MeV/sec) 

18 

.02-  .05 

MeV 

Gamma 

R, 

(MeV/sec) 

19 

Total  Gamma  Number 

Rate 

(Gaitunas/sec) 

20 

Total  Gamma  Energy 

Rate 

(MeV/sec) 
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This  report  deals  mainly  with  the  analysis  of  the 
235 

fission  of  U  (fission  spectra  fission) .  The  other  4 
fuels  covered  by  the  study  were  analyzed  and  the  results 
were  generally  the  same  as  those  for  U  (see  Table  II). 


Table  II 

Fuels  Covered  by  This  Study 
235 

U  (Fission  Spectra) 

238 

U  (Fission  Spectra) 

239 

Pu  (Fission  Spectra) 

235 

U  D  (14  MeV) 

U238  (14  MeV) 


Assumptions  and  Limitations 

A  large  amount  of  fission  yield  and  decay  data  was 
needed  for  this  study.  To  cover  all  the  possible  fission 
products,  one  would  have  to  include  data  on  approximately 
107  mass  chains.  This  would  involve  roughly  400  fission 
fragments  and  500-600  radioactive  fission  products.  Col¬ 
lection  and  compilation  of  data  on  all  possible  fission 
products  would  have  taken  more  time  than  was  allotted  for 
this  study.  The  data  used  was  limited  to  those  fission 
products  of  U  (fission  spectra  fission)  which  are  members 
of  mass  chains  with  a  fission  yield  of  1%,  or  higher.  This 
amounted  to  mass  chains,  which  involve  123  fission  fragments 
and  178  radioactive  fission  products.  Specifically,  data  for 
fission  products  in  mass  chains  84  through  105  and  129  through 
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149  were  collected  and  used  in  the  program.  As  a  result, 
the  program  accounted  for  the  following  percentages  of 
fission  products  from  the  respective  fuels  (see  Table  III) . 


Table  III 


Fission  Products 

Accounted  For 

Fuel 

Program  Product  Yield 

u235 

(Fission  Spectra) 

98.8% 

0238 

(Fission  Spectra) 

95.5% 

Pu239 

(Fission  Spectra) 

88.7% 

O235 

(14  MeV) 

82.9% 

u238 

(14  MeV) 

83.4% 

The 

program  will  compute 

results  for  any  specified 

amount  and  mixture  of  these  fuels.  A  fraction  of  fission 
products  accounted  for  from  the  specified  device  is 
computed.  The  rate  results  are  adjusted  by  multiplication 
by  the  inverse  of  this  fraction  so  they  will  reflect  an 
"adjusted"  100%  fission  product  coverage.  This  assumes  the 
gamma  radiation  from  the  missing  products  is  the  same  as 
that  from  the  included  products. 

The  most  serious  limitation  in  this  study  is  a  result 
of  the  fact  that  gamma  radiation  data  was  not  available  for 
many  of  the  short  lived  products.  Of  the  178  products  con¬ 
sidered,  gamma  radiation  data  was  not  available  for  84 
products.  Of  these  84  products,  19  had  half  lives  over  1 
minute,  6  had  half  lives  over  10  minutes,  and  2  had  half 
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lives  over  1  hour.  Using  the  fact  that  relatively  all 
activity  dies  out  after  10  half  lives  as  a  basis,  it  is 
considered  that  1000  sec  after  fission  is  a  conservative 
minimum  for  valid  program  energy  rate  results.  Thus, 
program  energy  rate  results  are  not  valid  for  the  time 
period  of  most  interest,  .1  sec  to  2  min  after  fission, 
however,  an  approximation  developed  in  Chapter  IV  allows 
total  gamma  energy  rate  calculations  for  this  time  region. 

Relative  gamma  intensities  were  available  for  16  of 
the  84  radioactive  products  with  unknown  absolute  gamma 
intensities. 

Table  IV  lists  gamma  radiation  data  for  2  fission 
products  to  illustrate  what  is  meant  by  absolute  and  rel- 


ative  gamma 

intensities 

(Ref 

8:224)  . 

Table  IV 

Kr89 

— 

— 

Rb90 

Gamma 

(MeV) 

Relative 

Intensity 

Gamma 

(MeV) 

Absolute 

Intensity 

0.23 

85.0 

0.53 

4% 

0.36 

28.0 

0.72 

4% 

0.51 

42.0 

0.83 

56% 

0.60 

100.0 

0.86 

6% 

others 

<  100.0 

others 

<  56% 

8  9 

For  the  decay  of  Kr  ,  only  relative  gamma  intensities  are 
known,  i.e.  the  0.60  MeV  gamma  is  the  most  frequent  gamma, 
and  for  every  100  0.60  MeV  gammas  there  would  be  42  0.51 


MeV  gammas.  The  number  of  any  gamma  produced  per  decay  of 
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89  90 

Kr  is  unknown.  For  the  decay  of  Rb  ,  absolute  gamma 

90 

intensities  are  known,  for  100  Rb  decays,  56  0.83  MeV 
gammas  would  be  produced.  It  was  observed  that  the  absolute 
intensities  of  the  most  frequent  gammas  from  all  the  products 
considered  generally  varied  from  100%  to  1.0%,  with  a  few 
exceptions  lower  than  1.0%.  Based  on  this  range  of  absolute 
intensity  values,  the  absolute  intensities  of  the  most  fre¬ 
quent  gammas,  from  the  16  products  for  which  only  relative 
intensities  were  known,  were  sat  at  1.0%.  With  the  absolute 
intensity  of  the  most  frequent  gamma  known,  absolute  inten¬ 
sities  of  the  remaii  ig  gammas  are  easily  calculated.  The 
absolute  intensity  of  1.0%  was  considered  to  be  generally 
low,  and  was  chosen  so  that  the  chances  that  estimated  gamma 
radiation  would  not  dominate  over  known  gamma  radiation  in 
the  results  would  be  minimized.  The  gamma  decay  data  for 
the  remaining  products  with  unknown  absolute  gamma  decay 
data  were  set  at  zero.  This  dictates  that  the  short  time 
(less  than  1000  sec  after  fission)  energy  rates  and  gamma 
number  rates  will  be  low. 

The  program  activity  results  (Rate  1,  see  Table  I)  do 
not  have  the  limitation  previously  described.  The  half  lives 
of  all  products  used  in  the  program  were  known.  However,  most 
mass  chains  have  some  members  for  which  half  lives  are  un¬ 
known.  This  study  used  the  first  member  of  the  mass  chain 
with  a  known  half  life  as  the  first  member  of  that  chain. 

Table  V  shows  an  example  of  the  mass  chain  information  used 
to  compile  data  for  the  program  (Ref  9:24). 
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Table  V 


Typical  Mass  Chain 


Nuclide 

Half  Life 

Cumulative  Fission  Yield 
(for  Fission  Spectra  U^5) 

34  Se96 

Unknown 

Not  Listed 

35  Br96 

Unknown 

7.56  x  10'9% 

36  Kr96 

Unknown 

9.08  x  10_5% 

37  Rb96 

.23  sec 

.0519% 

38  Sr96 

4.0  sec 

1.69% 

39  Y96 

2.3  min 

5.48% 

40  Zr96 

Stable 

6.22% 

96 

Thus,  37  Rb  was  used  as  the  first  member  of  the  chain 

with  a  fission  yield  of  .0519%  for  fission  spectra  U235. 

96 

The  second  member  was  38  Sr  ,  with  an  individual  yield 
of  1,64%  (1.69%  -  .0519%),  etc.  The  first  members  of  the 
mass  chains  with  known  half  lives  generally  had  half  lives 
of  .1  sec  to  1  min.  It  is  probable  that  the  members  with 
unknown  half  lives  have  very  short  half  lives  (less  than  1 
sec)  since  they  are  farther  from  stability  than  the  members 
with  known  half  lives.  Thus  some  activity  before  1  sec 
after  fission  is  not  accounted  for  and  therefore  program 
activity  results  are  considered  low  for  times  less  than  1 
sec  after  fission. 

Data  References 

The  references  used  to  gather  fission  yield  and  decay 
data  are  critical  to  the  validity  of  this  study.  They  are 
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listed  here  in  the  text  to  give  the  reader  an  idea  of  the 
scope  of  the  basic  data  used  in  this  study. 

Fission  yield  data  was  gathered  exclusively  from  the 
1972  General  Electric  Company  report,  "Compilation  of 
Fission  Product  Yields"  by  Meek  and  Rider  (Ref  7).  This 
report  presents  a  consolidated,  current  listing  of  this 
data. 

Half  life  and  decay  branching  data  was  gathered 
mainly  from  the  same  report,  with  occasional  reference 
to  other  sources  for  confirmation  purposes. 

Gamma  radiation  data  was  gathered  from  four  sources; 
"Handbook  of  Chemistry  and  Physics"  (Ref  10) ,  "Table  of 
Isotopes"  (Ref  8) ,  "Gamma-Ray  Energies  and  Intensities" 

(Ref  6) ,  and  applicable  editions  of  the  periodical.  Nuclear 
Data. 
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II.  Calculation  Method 

The  basic  inputs  usk d  in  this  study  are  device  yields 
device  composition  and  time  after  fission.  The  outputs 
are  various  energy  and  number  rates  as  described  in  Table  I. 

First,  the  number  of  atoms  of  each  fission  fragment 
produced  by  the  device  is  calculated.  The  device  yield  and 
device  composition  figures  are  used  to  compute  the  number 
of  fissions  of  each  type  fuel  that  occurred.  The  number 
of  fissions  of  each  fuel  is  then  multiplied  by  appropriate 
fission  yield  data  to  compute  the  number  of  atoms  of  each 
fission  fragment  produced  by  that  fuel.  The  total  number 
of  atoms  of  each  fission  fragment  produced  by  the  device 
is  then  calculated  by  adding  up  the  amounts  of  that  fission 
fragment  produced  by  the  component  fuels. 

Next,  the  initial  number  of  atoms  of  particular  fission 
fragment  and  the  specified  time  after  fission  are  used  to 
calculate  the  components  of  the  various  rates  described  in 
Table  I  caused  by  the  decay  of  this  particular  fission  frag¬ 
ment  and  its  daughters.  A  two  step  procedure  is  used  to 
compute  the  various  rates  for  each  fission  fragment  decay 
chain.  Running  totals  of  these  rates  are  kept  and  will 
reflect  output  rates  for  the  entire  device  after  the  pro¬ 
cedure  has  been  applied  to  all  fission  fragments. 

The  first  step  in  the  procedure  computes  the  number 
of  atoms  of  the  fission  fragment  and  its  daughters  present 
at  the  specified  time  after  fission.  Although  the  various 
isotopes  here  may  be  the  same  as  those  present  in  the 
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decay  chains  of  other  fission  fragments,  this  step  only 
accounts  for  those  particular  atoms  present  as  a  direct 
result  of  the  initial  quantity  of  the  particular  fission 
fragment  being  processed.  Thus  after  the  initial  quantities 
of  all  fission  fragments  have  been  processed,  rates  reflect¬ 
ing  the  effects  of  the  entire  device  will  have  been  calcu¬ 
lated.  Th3  basic  radioactive  transformation  equations, 
"Bateman  Equations"  (Ref  7:243)  were  used  in  this  step. 


Given  a  simple  decay  chain  for  fission  fragment  A^: 


A1  -►  A2  A3  A4  +  A5  +  etc. 


Aa(t)  -  A1(0)  e 


-Ait 


A2(t)  =Al(0)(A7A7e  +  Xrt  6 


Or  Generally: 


-A2fc 


-At  -A,t 

An(t)  =  Cl  e  +  c2  e  +  •  •  •  • 


Where:  = 


X^X~....X  i 
±  jl  n-i 


l2~a1'  '  3  •  •  *  •  A1 


AjtO) 


A.  X-.  .  •  *  A  « 

.’i  _  12  n-1 

l-2  a2-A!>  (X3-X1J 


VV  ai(o) 


Where:  t  is  the  specified  time  after  fission. 

A^(0)  is  the  initial  number  of  atoms  of  fission 
fragment  A^. 

A  (t)  is  the  number  of  atoms  of  the  nth  member  of 
n 

the  chain  present  at  time  t. 
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An  is  the  decay  constant  of  the  nth  member  of  the 
chain. 

These  equations  are  easily  adapted  to  computer  pro¬ 
cessing  with  proper  indexing  of  the  particular  variables. 

Decay  chain  branching  occurs  in  many  of  the  fission 
fragment  decay  chains.  For  the  simple  branching  decay 


the  decay  is  split  into  indexed  chains  as  indicated.  The 
branching  fractions  are  readily  available  in  the  references; 


i.e.  the  fraction  of  A2  that  decays  to  A^  3  (BF^)  and  the 
fraction  of  A2  that  decays  to  A2  3  (BF2> . 

The  following  equations  account  for  all  atoms  present 


at  time  t  for  the  simple  branching  decay  chain  shown  above. 


Ax(t)  =  Ax(0)  e 


-xifc 


(  X1  ~Xl^  X1  -X2^  \ 

Aa(t)  =Ai(0)  (17x7 e  +  7x^*7 e  )  (8) 


A,  „(t)  =  A, (0)BF 


(' 
«.  I  TT 


A  A  "^1^ 
A1A2  e  1 
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-Xjt 


A2f3(t)  =  Al(0)BF2 


X1X2 


_X2t 


X1X2 


~Xjt 


+  (X1-x2) (X3-X2)  +  (X1-X3) (X2-X3) 


(10) 


These  equations  are,  again,  easily  adapted  to  computer 
routines  with  proper  indexing  of  the  particular  variables, 
even  in  quite  complicated  ;  *anching  situations. 

The  second  step  in  the  procedure  multiplies  the  number 
of  atoms  of  each  radioactive  isotope  by  its  respective  decay 
constant  to  yield  activity  (decay/sec)  of  that  isotope  in 
the  chain.  This  activity  is  then  multiplied  by  decay  factors 
for  the  particular  isotope  to  give  the  rates  contributed  by 
this  isotope,  eg  activity  times  total  gamma  energy  per 
decay  of  this  isotope  would  yield  total  gamma  energy  per 
second  contributed  by  the  number  of  atoms  of  this  particular 
isotope  present  in  the  chain  oeing  processed.  The  rates 
from  all  isotopes  in  the  chain  are  then  summed  to  reflect 
the  rates  produced  by  the  entire  chain. 

Appendix  A  gives  a  listing  of  the  Computer  Program  used 
to  perform  the  calculations. 
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III .  Program  Results 

The  results  produced  by  the  program  consist  of 

various  energy  and  number  rates  as  described  in  Table  I. 

These  rates  are  quantitative  and  represent  the  output  from 

the  products  of  a  device  of  specified  yield  and  composition. 

The  results  presented  in  this  chapter  are  scaled  to  units  of 

235 

rate  per  fission  of  U  (fission  spectra  fission) .  These 

units  compare  more  readily  to  other  studies  and  eliminate 

the  need  to  specify  the  device  yield.  As  previously 

mentioned  in  Chapter  I ,  the  results  for  the  other  fuels 

235 

are  generally  the  same  as  those  computed  for  U 
Total  Gamma  Energy  Rate 

Figure  1  shows  the  total  gamma  energy  rate  as  a 
function  of  time  after  fission  as  calculated  by  the  computer 
program.  The  energy  rate  values  are  considered  low  for  time 
of  less  than  1000  sec  because  of  the  unknown  gamma  decay 
factors  on  many  short-lived  products  which  were  arbitrarily 
set  at  zero  in  the  program  data. 

Spectrum 

The  computer  code  produces  17  gamma-e~ orgy-group 
energy  rates  as  shown  in  Table  I.  These  results  are  subject 
to  the  same  limitations  as  the  total  gamma  energy  rate.  In 
Figure  2,  the  group  er  >rgy  rates  are  consolidated  into  3 
energy  groups  and  normalized  so  that  total  energy  rate 
equals  one.  Figure  2  indicates  very  erratic  spectral  be¬ 
havior  over  the  times  covered,  however  a  general  "softening" 
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TIME  AFTER  FISSION  (SEC) 

Figure  2.  Normalized  Gamma-Energy-Group  Energy  Rates 
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trend  (the  tendency  for  the  total  gamma  energy  to  be  composed 
of  mostly  lower  energy  gammas)  is  noted  as  times  progress 
beyond  1000  sec. 

Perhaps  a  better  look  at  the  gamma  radiation  spectrum 
is  given  by  Table  VI  which  lists  average  energy  per  gamma 
as  a  function  of  time  after  fission.  The  limitations 
below  1000  sec  still  a^-ply. 


Figure  3  shows  the  total  product  activity  as  a  function 
of  time  after  fission.  These  results  are  considered  accurate 
for  times  greater  than  1  sec  after  fission.  The  activity 
results  for  times  below  1  sec  are  considered  low  because 
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products  with  unknown  half  lives  (very  short  lived  products) 
were  not  included  in  the  computer  code. 
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IV.  Early  Time  Approximations 

Although  the  program  results  presented  in  Chapter  III 

are  considered  good  for  late  times,  these  results  do  not 

cover  the  time  period  of  most  interest;  .1  sec  to  2  min 

after  fission.  Two  methods  were  used  in  an  attempt  to 

find  a  procedure  that  would  produce  good  total  gamma  energy 

rate  values  for  the  .1  sec  to  2  min  time  period;  use  of  the 
-1  2 

familiar  "T  *  Law"  and  an  independent  approach  developed 

by  the  author.  As  in  Chapter  III,  the  results  presented  in 

235 

this  chapter  are  in  units  of  rate  per  fission  of  U  (fission 
spectra  fission) . 

-1  2 

"T  "  Law"  Approximation 

-1  2” 

The  familiar  "T  *  type  approximation  is  cited  as 

a  method  of  obtaining  product  radiation  rate  values  for 

times  less  than  1000  hours  after  fission  (Ref  4:421)  or  6 

months  alter  fission  (Ref  1:323).  It  is  stated  that  this 

approximation  is  good  only  for  times  in  excess  of  10  sec 

after  fission  (Ref  11:  10.H-1029).  This  method  involves 

the  use  of  the  equation: 

A(t)  =  Aot_1'2  (11) 

where:  A  is  the  radiation  rate  at  unit  time 
o 

after  fission. 

t  is  the  time  after  fission  of  interest 
in  units  compatible  with  Aq. 

A (t)  is  the  approximated  radiation  rate 
at  time  t. 

Thus.,  if  the  radiation  rate  at  any  time  after  fission  i3 
known,  an  approximation  of  the  rate  for  any  other  time  can 


20 


GNC/PH/74-3 


be  made  (within  the  specified  time  .limits) . 

Using  the  program  total  gamma  energy  rate  at  1000  sec 
after  fission  as  Aq,  this  approximation  was  used  to  compute 
rate  values  for  the  other  times  after  fission.  Figure  4 
shows  the  approximation  results  compared  with  the  program 
results. 

Both  the  approximation  results  and  the  program  results 
were  integrated  by  trapezoidal  numerical  methods  to  deter¬ 
mine  total  product  gamma  energy  released  per  fission.  The 

time  span  covered  by  the  integration  was  from  10  ^  sec 
12 

through  10  sec  after  fission.  The  integration  time  period 
starts  after  nuclear  weapon  burn  times  (considered  less  than 
10~7  sec) .  Negligible  gamma  energy  per  fission  is  contrib¬ 
uted  by  times  greater  than  105  sec  in  both  cases.  The  results 
were  308  MeV/fission  for  the  approximation  and  5.03  MeV/fission 
for  the  program  results.  It  appears  the  approximation  allows 
far  too  much  product  gamma  radiation  when  one  compares  the 
integration  results  to  generally  accepted  values  of  around 
6  MeV/Fission  (Ref  4:13).  Even  if  only  the  time  span  of 
interest,  .1  sec  through  2  min,  is  used  as  the  integration 

period,  the  approximation  results  yield  a  product  gamma 

••1  2 

energy  of  more  than  40  MeV/fission.  Clearly  the  "T 

Law”  approximation  produces  total  gamma  energy  rates  that 

are  too  large  over  the  .1  sec  to  2  min  time  after  fission 

-1  2 

region.  This  would  be  expected,  since  the  use  of  the  "T 
Law"  at  times  below  10  sec  is  not  considered  valid  (Ref  10: 
1029)  and  because  numerous  experimental  and  theoretical  total 
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gamma  rate  curves  show  "flattening"  tendencies  at  times 
less  than  1  to  100  sec  after  fission  (Ref  11:1328,  2:2-19, 
3:B811,  5:B822) . 

Independent  Approximation 

The  various  program  result  rates  were  analyzed  in  an 
attempt  to  find  trends  that  could  extend  the  long  time 
total  gamma  energy  rates  into  the  short  time  region.  The 
program  total  gamma  energy  rates  and  total  activity  rates 
were  used  to  compute  average  gamma  energy  per  decay  values. 
A  straight  line  approximation  fits  the  average  gamma  energy 
per  decay  calues  very  well  at  times  of  1000  sec  to  10,000 
sec  after  fission  as  shown  in  Figure  5.  This  time  regime 
was  selected  because  1000  sec  is  the  minimum  time  at  which 
the  program  total  energy  rate  results  are  considered  valid. 
The  straight  line  approximation  for  average  gamma  energy 
per  decay  is  given  by  the  following  equation: 

R(t)  =  Rot_‘104  (12) 

where:  Rq  is  the  average  gamme  energy  per  decay 
at  unit  time  after  fission. 

t  is  the  time  after  fission  of  interest  in 
units  compatible  with  Rq. 

R (t)  is  the  approximation  of  average 
gamma  energy  per  decay  for  time  t. 

Figure  5  shows  the  program  average  gamma  energy  per 
decay  values  and  approximation  values  based  on  the  program 
value  of  1.44  MeV/decay  at  1000  sec.  As  would  be  expected, 
the  approximated  values  are  larger  than  program  values  as 
time  goes  below  1000  sec  since  program  total  gamma  energy 
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rates  do  not  include  gamma  energy  from  some  short  lived 
product  decays.  As  previously  stated,  program  total  activity 
rates  are  considered  good  for  times  down  to  1  sec. 

The  approximated  average  gamma  energy  per  decay  values 
(MeV/decay)  were  then  multiplied  by  program  total  activity 
rate  values  (decays/sec)  to  produce  approximated  total  gamma 
energy  rate  values  (MeV/sec)  for  times  under  1000  sec.  The3e 
results  are  compared  to  program  results  in  Figure  6. 

It  should  be  noted  that  at  times  below  1  sec,  the  times 
when  program  total  activity  values  are  considered  low,  the 
approximated  average  gamma  energy  per  decay  values  are  in¬ 
creasing  above  2.5  MeV/decay.  Although  actual  cases  where 
more  than  2.5  MeV  of  gamma  energy  is  released  in  a  decay 
can  be  cited,  such  decays  are  rare.  Thus  while  average 
gamma  energy  per  decay  values  are  probably  high,  for  times 
less  than  1  sec,  and  the  total  activity  rate  is  considered 
low,  these  errors  will  tend  to  compensate  each  other,  to 
some  extent. 

In  Figure  7,  some  experimental  values  for  total  gamma 
235 

energy  rates  from  U  as  measured  by  Fisher  and  Engle  in 
1964  (Ref  2sB811)  are  compared  with  the  independent  approxi¬ 
mation  results.  It  should  be  noted  that  the  5  experimental 
values  fall  within,  and  pretty  much  cover,  the  .1  sec  to 
2  min  after  fission  period  of  interest  in  this  study.  The 
experimental  values  fit  the  approximation  results  very  well. 

Integration  of  the  approximation  results  over  the  time 
period  10~  through  10A  sec  results  in  a  total  gamma  energy 


25 


GAMMA  RATE  (MEV/SEC-FISSI0N) 


GNE/PH/74-3 


TOTAL  GAMMA  RATE  ( MEV/SEC-FISSION ) 


GNE/PH/74-3 
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Figure  7.  Total  Gamma  Energy  Rate  per  Fission 
Independent  Approximation  Compared 
with  Experimental  Values. 
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per  fission  value  of  8,75  MeV/fission.  This  is  close  to  the 
generally  accepted  value  of  6  ±  1  MeV  per  fission  (Ref  4s 13). 

Independent  Approximation  Compared  with  "Firefly" 

Approx ima t ion 

The  Air  Force  Weapons  Lab,  Kirtland  AFB,  N.  M. ,  is 
developing  a  code,  "Firefly",  to  calculate  weapon  product 
radiation  dose  to  aircraft.  A  memo  (Ref  12)  was  supplied 
to  the  author  describing  the  product  gamma  dose  calculations 
proposed  for  use  in  the  code. 

The  AFWL  approximation  uses  a  modification  of  the 
following  equation  (Ref  2:2-6)  as  a  source  of  energy  rate 
values. 

A(t)  =  c(l  +  t)-1'2  (13) 

where:  C  =  1.9 

t  is  the  time  after  fission  in  seconds., 

A ( t)  is  the  total  gamma  energy  rate 
(MeV/sec-fission)  at  time  t. 

This  equation  yields  a  total  gamma  energy  per  fission  value 
of  9.5  MeV  when  integrated  over  the  period  t  =  0  to  ®.  For 
use  in  the  Firefly  code,  equation  13  was  scaled  so  that  its 
integral  would  equal  7.7  MeV,  the  arithmetic  mean  of  9.5 
MeV  and  the  generally  accepted  value  of  6.0  MeV  (Ref  4:13). 
Thus  equation  3,  with  C  equal  to  1.55,  was  used  to  compute 
total  gamma  energy  rate  values  for  the  Firefly  code.  A 
comparison  of  the  results  of  Firefly  total  gamma  energy 
rate  values  with  the  author ' s  independent  approximation 
results  is  shown  in  Figure  8. 

To  evaluate  the  effects  these  two  approximations  would 
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Figure  8.  Total  Gamma  Energy  Rate  per  Fission  - 

Independent  Approximation  Compared  with 
"Firefly"  Approximation. 
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have  on  dose  calculations,  both  were  integrated  over  in¬ 
creasing  time  periods.  Accumulated  total  gamma  energy  vs. 
time  was  plotted  for  both  methods  (see  Figure  9) .  The 
independent  approximation  integration  values  differ  by  as 
much  as  29%  at  20  sec  from  the  "Firefly"  integration 
values. 

The  Firefly  calculations  use  a  fixed  gamma  spectrum 
based  on  the  following  equation  (Ref  12:1): 


dN (E)  _  ,  „  _-l • IE 
~di —  1-1  e 


(14) 


where:  N (E)  is  the  number  of  gammas  with  energy  E. 

E  is  the  gamma  energy  in  MeV. 

Although  Firefly  calculations  used  different  gamma  energy 
groups  than  this  study,  integration  of  equation  14  over 
appropriate  energy  ranges  easily  illustrates  the  Firefly 
spectrum  in  terms  of  the  17  energy  groups  used  in  this 
study  (sec  Table  I) .  Figure  10  compares  the  Firefly  fixed 
spectrum  with  the  spectrum  calculated  in  this  study  at  1000 
sec  after  fission.  The  gamma  number  values  used  are  scaled 
so  that  total  number  of  gammas  over  all  groups  equals  .978, 
the  integral  of  equation  4  over  the  energy  range  .02  to  80 
MeV  (the  range  of  energy  values  covered  in  the  17  groups  used 
in  this  study) .  It  should  be  noted  that  the  average  gamma 
energy  over  this  range  of  energies  is  .848  MeV  for  the  "Fire¬ 
fly"  spectrum  and  1.054  MeV  for  the  results  of  this  study  at 
1000  sec  after  fission. 
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•’iyure  10.  "Firefly"  Spectrum  Comparison 
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V.  Conclusions 
235 

For  a  device  using  only  U  (fission  spectra)  as  a 

fuel,  some  of  the  results  of  this  study  can  be  used  with 
confidence.  Total  activity  values  can  be  used  for  times 
after  fission  greater  than  1  sec.  Gamma  group  energy  rates 
can  be  used  for  times  greater  than  1000  sec  after  fission. 
The  program  results  can  be  used  for  total  gamma  energy 
rates  for  times  down  to  1000  sec  after  fission.  The 
independent  approximation  results  can  be  used  for  total 
gamma  energy  rates  at  times  less  than  1000  sec  after  fission 

Although,  as  previously  mentioned,  the  results  for 

235 

devices  composed  of  fuels  other  than  U  (fission  spectra) 

were  generally  the  same,  it  is  felt  the  lower  fraction  of 

products  included  in  the  program  for  these  fuels  (see  Table 

III)  makes  their  results  subject  to  more  question  than  for 
235 

the  case  of  U  .  The  results  for  other  fuels  might  be 
slightly  different  if  the  missing  fission  products  were 
accounted  for  in  the  program  data. 

If  complete  data  could  be  obtained  for  the  short 
lived  fission  products  and  added  to  the  program  data,  the 
program  could  produce  very  accurate  results  for  short 
times  after  fission.  The  calculation  methods  used  in  this 
study  are  simple  and  analytical.  The  results  produced  by 
this  method  are  only  limited  by  the  accuracy  and  complete¬ 
ness  of  the  fission  yield  and  decay  data  used  in  the 
calculations. 
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Appendix  A 


The  computer  program  used  in  this  study  was  written 
in  FORTRAN  IV  language  for  use  on  the  AFIT  CDC-6600 
computer.  The  main  program  (GAMA)  and  the  subroutine 
(BATMN)  listings  are  included  in  this  report.  The  two 
data  files  containing  fission  yield  and  decay  data  are 
not  included*  They  are  considered  too  lengthy  and  to  be 
of  little  interest  to  the  majority  of  those  who  will  read 
this  report*  A  copy  of  the  entire  code  is  on  file  at  the 
Physics  Department  (AFiT-ENP) ,  Wright-Pat  ter  son  AFB,  Ohio. 
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******************  PROGRAM  GAMA  PAGE  1  ***************** 


p*  ***************************  *•****•■  ******  *******  *********** 

C*  * 

C*  PROG® AM  GAMMA  COMPETES  ANC  PRINTS  CUT  VARIOUS  ENERGY  * 
C*  AND  NUMRER  RATFS  AS  QPSCRIBEO  IN  TABLE  I  FOR  A  FISSION  * 
C*  OEVICE  OF  SPECIFIC  YIFIC  ANC  COMPOSITION.  THE  PROGRAM* 
C*  13  DIMENSIONED  SO  THAT  CITFUT  FOR  IP  TO  50  TIMES  AFTER  * 
C*  FI SSI OM  HILL  PE  PRODUCED.  * 

C*  * 

C*  TAPE  5  MUST  BE  ATTACHED.  THIS  FILE  IS  COMPOSED  OF  * 
C*  OECAY  DATA  (NC,  N,  NO*  NFC.  N80,  FF,  GG,  ANO  PC)  FOR  * 
C*  FISSION  FRAGMcNTS.  THE  FILE  U3£o  FOR  THIS  STUDY  * 

C*  INCLUDED  DATA  FOR  12-*  FISSICN  FRAGMENTS.  ALL  FISSION  * 
C*  FRAGMENTS  ARE  INDEXED  ANC  FISSION  FRAGMENT  OECAY  OATA  * 
C*  MtJST  BE  ARRANG'D  IN  NUMERICAL  ORDFR  ACCORDING  TO  THE  * 
C*  INDEX  NUMBERS.  * 

C*  * 

C*  PROGRAM  CONTROL  PARAMETERS  (MIME,  TA  ,  YIELD,  ANO  YCOMP)  * 
C*  ARE  ATTACHED  AS  A  OATA  DECK  ANO  READ  IN  FROM  THE  FILE  * 
C*  "INFUT“.  * 

C*  * 

C*  SUBROUTINE  BLOCK  DATA  MUST  PE  ATTACHEO.  IT  SETS  VALUES* 
C*  FOR  (YFF,  NIO.  ANO  FPKT).  * 

C*  * 

C*  SUBROUTINE  BATMN  MUST  BE  ATTACHED.  IT  CALCULATES  THE  * 
C*  RATE  CONTRIBUTIONS  OF  EACH  FISSION  FRAGMENT  AS  * 

C*  DESCRIBED  IN  CHAPTER  II. 

C*  * 

C*  PROGRAM  EXECUTION  REQUIRES  37K  ANO  120  SEC  TO  COMPUTE  * 
C*  OUTPUT  RATES  FCR  20  TIMES  AFTER  FISSICN.  * 

C*  * 

PROGRAM  GAMft( INPUT, OUTPUT, TAPE5, PLOT) 

DOUBLE  PRECISION  GGTT (20, E 0) , GGT (2 0) 

DIMENSION  TA (52) .YCOMP  (7) ,FFKT (7 ) , YFF ( 123 , 7) , FF 0 ( 123) , 
2NIO (123 ) ,GGTF(52> ,N<3) , H 9 1  3) ,NFC <3 ) ,NPO ( 3) ,NBP(3,3> , 
3EF( 3,3) ,00(10,3) ,GG (20,?, 3> 

COMMON  FF01,T,NC,N,Nn,NFC, NBO, N8F, 3F,OC,GG,GGT 
COM  MON/ FFQ/  TF, NIO, FPKT 
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******************  FROGRAM  oAMA  CAGE  2  ***************** 


10  FORMAT  (15) 

11  FORMAT  (3E1".3) 

20  FORMAT  (8E10.4) 

70  FORMAT  (*1*, ’TOTAL  GAmkm  PATE*) 

71  FORMAT  (’-’,’TIME  AFTfF  0LAST*,5X, ’GAMMAS  R'R  SEC*, 

213X,*MEV  PER  SEC* , 1£X  , ’ACTIVITY  (DECAYS  PEP  SEC)*) 

72  FORMAT  (*  *) 

74  FORMAT  (1X,1PE14.7,7X,1PC23.15,5X,1°C23.15,5X,1RD23.15> 

80  FORMAT  < *-* ,5 3< *-*) , ’FISSION  FRAGMENT  YIELDS*, 53 (*-*) ) 

31  FORMAT  (*-*,*TOTAL  MJMeEP  OF  FISSIONS  * , 1PE 14. 7 ,5X, 
2*RATIO=  ** 1FE14.7) 

82  FORMAT  (*-*,*FF  ID* , 10X,* YIEL D  OF  FF  (NUMPfcR )* ,5 X, 

2*Y IELO  OF  FF  (FRACTION  CF  TOTAL  FISSIONS)*) 

83  FORMAT  (*  *) 

85  FORMAT  (IX, l4,*-*,Ifc,*G’,5X, 1PL14.7, 11X,F13. 8) 

87  FORMAT  ( IX, 14 , ♦-*, 16, *M* , SX, l^E 14. 7, 1 IX ,F13. 8) 

30  FORMAT  ( *1* ,57 (*-*) , ’HEAFON  OAT A* , 57 ( *-*) ) 

91  FORMAT  (*-* ,*TOTAL  YIELO  *,1PE10.3,*  KILO  TONS”) 

32  FORMAT  (*-*,»FUEL  FRACTIONS*) 

93  FORMAT  (*-♦ ,*U235 (FS ) * ,7X ,*U236 (F S) * , 7X ,*PU2 33 (FS) *, 
26X,*U?35114M£V)*,4X,*U<38 (14MEV)* , 4X, ’PU239 ( 14  Me  V) *, 

34X , *  OTHER*, 3X) 

94  FORMAT  ( IX ,  7(  1PE10  .  ,3 ,SX) ) 

150  FORMAT  ( 8 1 1 J ) 

151  FORMAT  (8 F^.7) 

152  FORMAT  (4C20.13) 

200  FORMAT  («1C-".0U=  GAMMA  PATFS*) 

201  FORMAT  (*-TlM^  AFT  0LAST*,3X,*  *  ,  12X  ,*  G=>OUP  2*, 

212X,»G°0UP  3*  ,  12X, *C-RCLF  *■*»  12X  ,*  GR  )L'F  v*  ,  12X,  *G=-OUF  6*) 

233  FORMAT  (IX  ,  IP  "1  4. 7 , 24X  ,5  ( 1=31  4.  7,  ?K  )  ) 

204  FORMAT  (*-TIM_  AFTEF  ELAST*, 3Y , ’GROUP  7* , 12X ,’ 5RCUP  3*, 
212X, ’GROUP  M*,12X,*Gk0LP  10* , 11 X, *GROUP  11*, UX, 

3*GROUP  12*) 

206  FORMAT  (’-TIME  1FTEP  RLAST*,  3X  , *G ’OUP  1**,11X, 

2*G  RO'IP  14* ,  11 X ,  ’GROUP  lr->  ,11X  ,  *GR  1U°  16’,11X, 

3 ’GROUP  17* , 1 1 X , *G°OUP  18* ) 

208  FORMAT  ( IX , 1PE14.7  ,5X , t  ( 1F314.7 ,5X) ) 


****** 

* 

****** 

RE AO  10.NTIME 

READ  11, (TA (I) ,I=1,NTI“E> 

PE  AO  20.YIELO, (VCOMP(I),I=t,7) 


£**«•• »*  *«***•»*•*•******«•*  ********* **•**••******••*» 
C* 

C*  Rc AO  IN  PROGRAM  CONTROL  FARAMFT£R  S» 

C* 

c*  ********  *****.., 
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+« * **#♦♦♦**++*»++ *  PROGRAM  GAMA  PAGt  3  *♦*****♦***♦**♦♦♦ 


£4444* *4 44 ¥+4*44+ +4 4*44444444#***#+*+ *4 *44444+4*4+ *«*♦++♦«*+ 

c*  * 

C*  COMPUTE  THc  YIELD  (NUM8E*  OF  ATOMS)  OF  EACH  FISSION  ♦ 

C*  FRAGMENT  P^ODUCEC.  * 

C*  * 

Q4+4444+++44++4+++4++44++4++44+++*************************** 

00  25  J=l,123 
FFO(J)  =0, 

DO  2&  1=1,7 

25  FFO(J)  =FFO ( J) +YFF ( J, I) *YC0MP (I)  *Y IELD*FPKT (I) 


Q  +  |f,f  #4  +  ***  +  ♦*  *+>+  +  +  ♦♦*++**  *+**+++* 

C*  * 

C*  COMPUTE  THE  RATIO  OF  THE  ACTUAL  TOTAL  NUMBER  OF  FISSION* 
C*  FRAGMENTS  (2-  TIM£S  THE  TOTAL  NUMBER  OF  FISSIONS)  TO  THE* 
C*  COMPUTED  TCTAL  NUMBER  OF  FISSION  FRAGMENTS,  * 

C*  * 

Q4 444 + 4 *++ 4 +++*+++++++++++++++*+++++*++++*++*+*+++++++++ + 


FFOTOT-O. 

00  30  J=1 , 123 
30  FFOTOT=FFOTCT+FFO( J) 

FTOT  =  0. 

00  3  5  3=1,7 

35  FT OTsFTOT ♦YIELD* YCOMF(w)*FPKTCJ) 

PA  ~T0=  (2,*F  TOD/FFOTOT 


c*  * 

C*  PRINT  OUT  PAPAMErERS  DESCRIBING  THE  DEVICE#  * 

c*  * 

044*  4*** **♦♦*+♦ +*44**+4+*4**+++*+*+  4  44  ++  *4+  4 

PPINT  90 
PRINT  91, YIELO 
PRINT  9? 

PRINT  93 

PRINT  9A,YCC«P<1)  ,YCCM.F(2>  ,YC0MF(3) ,YC0MP(4> ,YC0MP(5> , 
2YC0MP(u) ,YCQHP(7) 

50  PRINT  61 , FTOT , DS  TIQ 
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#»*+«♦** *♦*♦***♦*  *  p^rGRAM  GAMA  PAGE  4  ♦  ♦♦♦**♦*♦♦♦♦***♦«• ' 


*****  *********************************  4*  ******** 

c*  * 

C*  A  LIST  INCLUDING  r ACM  FISSION  FRAGMENT  INOFY  NUMBER,  ♦ 

C*  FISSION  FRAGMENT  IDENTIFICATION  (ATOMIC  NUMRFp,  MASS  * 

C*  NUMPERt  AND  STATE),  DE  VICE  YIELO  CF  THAT  FRAGMENT  * 

C*  (NUMBER  OF  ATOMS)  AND  FRACTIONAL  YIELD  IS  PRINTED  OUT.  * 

C*  * 

Q*J ,******4************************************************** 

PRINT  80 
PRINT  82 


PRINT  83 


no 

84 

1  =  1. 

123 

PC-CFFO (I) 

/FTOT) *. 

1 

IF 

(I 

•  EO. 

33) 

GO 

TO 

86 

IF 

(I 

.EQ. 

37) 

GO 

TC 

86 

IF 

(I 

•  EQ.. 

40) 

GO 

TO 

86 

IF 

(I 

•  EC. 

46) 

GO 

TO 

66 

IF 

(I 

•  EO. 

?3) 

GO 

TO 

86 

IF 

(I 

•  EO. 

57) 

GO 

TO 

96 

IF 

(I 

.FG. 

90) 

GO 

TO 

86 

IF 

(I 

.EC. 

95) 

GO 

TO 

66 

IF 

(I 

•  EG. 

109) 

GO 

i  TC 

66 

IF 

(I 

.EC. 

112) 

GO 

TO 

6b 

IF 

(I 

•  EO. 

115) 

GO 

TO 

86 

FRINT 

85,1 

,NIO( 

I). 

FFO 

(I) 

GO 

TO 

84 

PRINT 

87,1 

,NID( 

I)  . 

FFO 

(I) 

84  CONTINUE 


c*  * 

C*  RFA9  FISSICN  FRAGmfNt  DECAY  PARAMETERS  FOR  A  FISSION  * 

C*  F 3 AG ME NT ,  COMMUTE  RATF  CCMTP IBUTIONS  OF  T«E  FRAGMENT  * 

L*  ANO  ADO  THEM  TO  RUNNING  TOTALS.  REPEAT  FOR  EACH  ♦ 

C ♦  FISSION  FRAGMENT  *  REPEAT  THE  PROCESS  FOR  EACH  tjmc  ♦ 

C*  AFTER  FISSION.  NOTE™ ALL  INITIAL  FISSION  FRAGMENT  * 

C*  NUMBERS  A°E  MULTIPLIED  PY  RATIO  SC  AN  ADJUSTED  1307.  * 

C*  FISSION  YIELD  COVERAGE  WILL  BE  REFLECTED  In  THE  OUTPUT. ♦ 

C*  ♦ 

Q***** ********************************************  ********** 
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♦  *  4  *  ♦*  *** *  *****  *  PROGRAM  GAMA  PAGE  5  ***************** 


DO  65  K=i , NTI ME 
TsTA(K) 

00  55  J=i,20 
55  GGT(J)=0. 

00  160  K*=  1 » 1 30 
FF01=FF0(KK)*RATIO 

READ  (5,150)  NC, (N(I)  ,1=1,3) , (NE(I) ,1  =  1, 3> 
IF  (EOF  (5)  ,N£.  0)  GO  TO  170 
PE AO  (5,150)  ( NCC ( I )  ,1  =  1,3), (N 30(1), 1  =  1, 3) 
00  140  JslfNC 

140  READ  (5,150)  (N8P(I ,J) ,1*1,3) 

00  141  J=l, NC 

141  READ  (5,151)  (BF(I,J)  ,1  =  1,3) 

00  1  42  J=  1,  NC 

NN=M(J) 

DO  142  JJ=1 ,NN 

142  RE  AO  (5,151)  (GG ( I , J J, J) , 1=1, 20 ) 

DO  143  J  =  1 ,  NC 

NN=N(J) 

143  READ  ( 5, 152 )  (OC (I , J )  ,  1=1 ,NN) 

CALL  3ATMN 

160  CONTINUE 
170  REWIND  5 

00  64  J  =  1 , 2  0 

64  GGTT ( J, K) =GGT ( J) 

65  CONTINUE 


Q++**»**+«*#*«+******«**#«+«#4###***««************ ********** 

0*  * 

C*  PRINT  THE  C0MPUT5C  RATES  FOP  EACH  TIME  AFTER  FISSION,  * 

C*  * 

£****+*»**+**«****»»********«******************************* 


73 


PRINT  70 
PRINT  71 
PRINT  72 
00  73  I=1,NTIME 

PRINT  74  ,TA  (I ) , GGTT  (19,1 ) ,GGTT (20fI) , GGTT (1,1 ) 
PRINT  200 
PRINT  201 
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»++*+#+»+»»+****++  PROGRAM  GAMA  PAG  £  6  *♦****♦****♦♦♦♦*♦ 


PRINT  72 

DO  202  1  =  1* NT  IMF 

202  PRINT  203,TA(I) * (GGTT  <  J,  I  > ,  J  =  2  ,  C) 
PRINT  200 
PRINT  204 
PRINT  ?2 

DO  205  1  =  1  *  NT  IMG 

205  PRINT  208, TA( I) , (GGTT( J,I) *J=7, 12) 
PRINT  200 
PRINT  206 
PRINT  72 

00  207  1=1, NT I ME 

2  07  PRINT  208 ,T A ( I ) , (GGTT ( J , I ) , J= 1 7 ,1 3) 
STOP 
END 
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Program  GAMA  Variable  Definitions: 


NTIME 

-  The  number  of  times  after  fission  for 
which  results  are  desired. 

TA  (I) 

-  The  times  after  fission  in  sec. 

Yield 

-  Total  fission  yield  of  the  device  in 
kilo-tons  TNT. 

YCOMP(I) 

-  The  fraction  of  the  total  yield  from 
the  Ith  fuel.  Table  II  lists  the 
fuels  in  order,  eg  fission  spectra 

235 

U  J  is  fuel  1. 

YFF(I,J) 

-  The  fission  yield  of  fragment  I  for 
the  fission  of  fuel  J. 

NID(I) 

-  The  identification  (atomic  number  and 
weight)  of  fission  fragment  I. 

FPKT ( J) 

-  The  number  of  fissions  of  fuel  J 
needed  for  one  kilo-ton  of  yield. 

GGTT (I , J) 

-  The  output,  ie  the  rate  I  for  time  J 
after  fission. 
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***♦¥*♦***♦**♦*♦♦♦  SUBROUTINE  0ATMN  PAGE  i  ************** 


£  +  *  +  +  +  * JMMMMMMMMMF****  +  <MMMHMMMMMMMHMM*'»  <MHHMMMMHMHMM* ****+*««* 


c*  * 

C*  SUBROUTINE  BATHN  COMPUTES  THF  RATE  CONTRIBUTIONS  OF  A  * 
C*  FISSION  FRAGMENT  AND  ADDS  THFM  TO  RUNNING  TOTALS  FOR  ♦ 
0*  THE  TIME  AFTER  FISSION  BEING  PROCESSED.  * 
C*  * 
C*  THE  INTIAL  NUMBER  OF  ATOMS  OF  THE  FRAGMENT,  THE  TIME  ♦ 
C*  AFTER  FISSICN,  AND  FRAGMENT  CECAY  PARAMETERS  ARE  THE  * 
C*  INRUT  S  TO  SUBROUTINE  BATMN.  * 
C*  * 


£»«»*** ♦**+***♦*****#***#+****#+♦♦+*#* ***4****************** 


SUBROUTINE  BATMN 

DOUBLE  PRECISION  GGT (20 ) , GGS <20 ) ,GGSS < 20 > 
DIMENSION  N(3> ,  NB  ( 3  >  , MFC <  3 )  ,DC  <  1  0, 3) , GG <20, 10  ,3) , 
2  N  RO ( 3 ) , NBP(3,3) ,  BF ( 3, 3) »FX ( 10) , C (1 0) 

COMMON  FFO,T,NC,N,NB,NFC,NEO,NBF,  BF, DC ,GG,GGT 

CO  1  1=1,20 
1  GGS  ( I )  =  0. 


«* ***»***«*+**« ***** 


c*  * 

c*  THE  INPUTS  ARE  USEC  TO  SET  UF  DECAY  EQUATIONS  (AS  * 
C*  ILLUSTRATED  BY  EQUATIONS  7-10)  FOR  THE  PARTICULAR  * 
C*  FRAGMENT  AND  ITS  DAUGHTERS.  THE  CALCULATIONS  ARE  THEN  * 
C*  MAQl  AS  DESCRIBE"  IN  CHAPTER  II  TO  COMPUTE  THE  RATE  * 
C*  CONTRIBUTIONS  FOR  THIS  FRAGMENT.  * 
C*  * 


00  200  11*1, NC 
N 1=  N  ( 1 1 ) 

N2=NR(II> 

N3=NRO(II) 

DO  4  I  =  1 , N1 

4  EX( I) =  EXP  <-DC( I f II) *T) 

DO  3  1=1,20 
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»+*»»+»»+*+**+****  SU3R0UT IK E  BATMN  PAGE  2 


?  GGSS<I>=0. 

00  3  K=N2, N 1 
00  9  1*1,  K 
Ca)=0C(K,II)  *FF 0 
00  17  M*1 , N 3 
N4=N0PCM, II) 

IF  (K  .GT.  N4)  C(I)=C(I)*FF(Hf II) 

IF  (K  .20.  N4>  C(I)=C(I)*FF<MfII) 

17  CONTINUE 

IF  (K  .EO.  1)  GO  TO  9 
00  5  J-l ,  K 

IF  (J  .EO.  K)  GO  TO  25 
IF  (J  .EO.  I)  GC  TO  20 

c<i)=c<i)*DC(j,ii>/  (occj,  in -oca  ,n>  > 

GO  TO  5 

20  C<I)=C<I)*DC<J,II) 

GO  TO  5 

25  IF  (I  .EO.  K)  GO  TO  5 

C(I)=C(I)/(UC(  J.H)  -DC  (I,  II)  ) 

5  CONTINUE 

9  CONTINUE 

A=0. 

DO  6  J  =  1 .  K 
6  A=  A+C  ( J)  #EX  ( J) 

00  7  J  =1  *  2  0 

7  GGSS<J)=GGSS{J)*GG  ( J  » K  ,  II )  *A 

6  CONTINUE 

00  121  1=1.20 

121  GGS<I)=GGS(I>  +GGSS  ( I ) 

200  CONTINUE 


C* 


* 


C+  THE  FRAGMENT  -ATE  CONI R I EUT ICNS  APE  ADDED  TO  PUNNING  ♦ 
C*"  TOTALS  FOP  THE  TIME  AFTEP  FISSION  3EING  PROCESSED.  * 

C*  * 


*#«+«++»»+* ********** 


DO  122  1=1,20 

122  GGT (I) =GGT ( I) +  GGS  (1) 
RETURN 
END 
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Subroutine  BATMN  Variable  Definitions: 

FFO  -  The  initial  number  of  atoms  of  the 

fission  fragment  created  by  the  device. 

T  -  The  specified  time  after  fission  (sec) . 

GGT(I)  -  A  running  total  of  rate  I:  eg 

GGT(I)  is  the  total  activity  rate 
(decays/sec) . 

GGS (I) ,GGSS (I)  -  Internal  subroutine  summations  of 

rate  I. 

NC  -  The  number  of  chains  in  the  decay  of 

this  fission  fragment  (See  Appendix  A) . 

N (I)  -  The  number  of  radioactive  members  along 

the  Ith  chain. 

NB(I)  -  The  first  member  of  chain  I  that 

branches. 

The  chain  that  chain  I  branches  off  of. 

The  number  of  branching  members  along 
chain  I. 

NBP(J,I)  -  The  Jth  member  of  chain  I  that  branches. 

BF(J,I)  -  The  branching  fraction  of  the  Jth 

member  of  chain  I  that  branches. 

DC(J,I)  -  The  decay  constant  of  the  Jth  member 

of  chain  I  (sec  j . 

-  The  Kth  decay  factor  of  the  Jth  member 
of  chain  I:  eg  GG  (20,  3,  2)  would  be 
the  total  gamma  energy  produced  by 
one  decay  of  the  3rd  member  of  the 
2nd  chain. 


NPC(I) 

NBO(I) 


GG(K,J,I) 
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